Zinc Impurities in the 2D Hubbard model 
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We study the two-dimensional Hubbard model with nonmagnetic Zn impurities modeled by bi- 
nary diagonal disorder using Quantum Monte Carlo within the Dynamical Cluster Approximation. 
With increasing Zn content we find a strong suppression of d-wave superconductivity concomitant 
to a reduction of antiferromagnetic spin fluctuations. Tc vanishes linearly with Zn impurity concen- 
tration. The spin susceptibility changes from pseudogap to Curie- Weiss like behavior indicating the 
existence of free magnetic moments in the Zn doped system. We interpret these results within the 
RVB picture. 
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Introduction Chemical substitution in high- Tc super- 
conductors provides a powerful probe into the complex 
nature of both the superconducting and normal state of 
these materials. Experiments substituting different im- 
purities for Cu show that nonmagnetic impurities (Zn, 
Al) are just as effective in suppressing superconductiv- 
ity as magnetic dopants (Ni, Fe) Based on Ander- 
son's theorem Q these results are taken as strong evi- 
dence for an unconventional pairing state described by 
an anisotropic order parameter with nodes on the Fermi 
surface. Indeed, recently the characteristic fourfold sym- 
metry of the (ia;2_y2-wave order parameter was observed 
in the spatial variation of the local density of states near 
the Zn impurity using scanning tunneling microscopy . 

Nuclear magnetic resonance experiments show that 
even a nonmagnetic impurity substituted for Cu induces 
an effective magnetic moment residing on the neighboring 
Cu sites [Q. In addition, the bulk susceptibility in im- 
purity substituted underdoped cuprate superconductors 
shows Curie- Weiss like behavior irrespective of the mag- 
netic structure of the dopants [Q. This indicates the ex- 
istence of free magnetic moments in the impurity doped 
Cu02 planes. The formation of these moments with sub- 
stitution of nonmagnetic impurities is explained to arise 
from breaking singlets |^ in an antiferromagnetically 
correlated host by removing Cu spins. 

In this letter we focus on the suppression of supercon- 
ductivity by Zn impurities and the change of the bulk 
magnetic susceptibility to a Curie- Weiss like behavior 
in underdoped systems. Despite a wide variety of the- 
oretical studies addressing these issues, a complete un- 
derstanding of the effects of Zn doping has not been 
achieved. These approaches are mostly based on the de- 
scription of impurities embedded in a BCS host |^ or 
on generalizations of the Abrikosov-Gorkov equations for 
nonmagnetic impurities in unconventional superconduc- 
tors using phenomenological pairing interactions |7|. In 
order to capture effects like moment formation by substi- 
tution of nonmagnetic impurities, it seems necessary to 
describe correlations and the scattering from impurities 
on the same footing within a microscopic approach. 

The most widely applied technique developed to de- 
scribe disordered systems is the Coherent Potential Ap- 



proximation (CPA) S. The CPA shares the same mi- 
croscopic definition |9[| as its equivalent for correlated 
clean systems, the Dynamical Mean Field Approxima- 
tion (DMFA) [0 . Both approaches map the lattice sys- 
tem onto an effective impurity problem embedded in a 
host that represents the remaining degrees of freedom. 
This single-site approximation neglects interference ef- 
fects of the scattering off different impurity sites (cross- 
ing diagrams) and correlation effects become purely local. 
Therefore, it inhibits a transition to a state described by 
a nonlocal (d-wave) order parameter and thus is not ap- 
propriate for our investigations here. 

The Dynamical Cluster Approximation (DCA) |^ 
[T^ , [T^ systematically incorporates nonlocal corrections 
to the CPA/DMFA by mapping the lattice system onto 
an embedded periodic cluster of size TVc- For Nc = 1 the 
DCA is equivalent to the CPA/DMFA and by increasing 
the cluster size Nc the length-scale of possible dynamical 
correlations can be gradually increased while the DCA so- 
lution remains in the thermodynamic limit. In the clean 
limit the DCA applied to the Hubbard model has been 
shown to describe the essential low energy physics of the 
cuprates jl^, |T^: It captures the antiferromagnetic 
phase near half filling and the transition to a supercon- 
ducting phase with dx2-y2-wa,ve order parameter at finite 
doping. In the normal state it exhibits non-Fermi liquid 
behavior in form of a pseudogap in the density of states 
and a suppression of spin excitations at low temperatures 
in the underdoped regime. 

In this letter we study the two-dimensional (2D) Hub- 
bard model including a potential scattering term accord- 
ing to the chemistry of Zn impurities using the DCA. 
We will show that potential scattering by Zn impuri- 
ties strongly suppresses superconductivity as well as spin 
fluctuations and changes the magnetic susceptibility to a 
Curie- Weiss like temperature dependence in the under- 
doped region. 

Formalism The correlated electrons in a Cu02 plane 
doped with Zn impurities are described by the 2D Hub- 
bard model with diagonal disorder 
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where we used standard notation. The disorder induced 
by the Zn sites occurs in the local orbital energies which 
are independent quenched random variables distributed 
according to some specified distribution P(ei,...,ejv) = 
Y[i=i P{^i) where N is the number of lattice sites. For a 
concentration x of Zn impurities we use the binary alloy 
distribution 

P(e,;) x5{e, + V/2) + (1 - x)6{u - V/2) , (2) 

where V is the energy difference between the Cu and Zn 
3d orbitals. For simplicity we use site- independent values 
for the hopping integral t and the Coulomb repulsion U. 

The microscopic derivation of the DCA algorithm was 
discussed in detail for correlated systems in |l^ and 
for disordered systems in j9|. The DCA has a simple 
physical interpretation for systems where the intersite 
correlations have only short spatial range. The corre- 
sponding self-energy may then be calculated on a coarse 
grid of Nc — selected K points only, where Lc is the 
linear dimension of the cluster of K points. Knowledge 
of the momentum dependence on a finer grid may be dis- 
carded to reduce the complexity of the problem. To this 
end the first Brillouin zone is divided into Nc cells of size 
{2tt/Lc)^ around the cluster momenta K. The propa- 
gators used to form the self-energy are coarse grained or 
averaged over the momenta K-|-k surrounding the cluster 
momentum K 

G{K,iu;n) = ^^G(K + k,zc^„) 

k 



N ^ iuJn - CK+k - S(K, iuJn) ' 

In Eq.(^) we make use of the fact that the approxima- 
tion of the lattice self-energy S(k) by the cluster self- 
energy S(K) optimizes the free energy S(K) 
summarizes the single-particle effects of the interaction 
term (second term) and the disorder term (third term) 
of the Hamiltonian, Eq.(|l]). The dispersion Cj^^j^ denotes 
the spectrum of the kinetic part (first term) of Eq. (|l|) . In 
order to avoid overcounting of self-energy diagrams on 
the cluster, the cluster excluded propagator 

g(K)-i =G(K)-i+E(K) (4) 

is used as bare propagator for the cluster problem. To 
diagonalize the disorder part it is convenient to per- 
form a Fourier transform to cluster real space Gij — 
^/NcJ2K^^^'^'^'~'^^^Si'K.) where the disorder is diago- 
nal. The inverse bare cluster propagator for a particular 
disorder configuration {ei} is then given by 

[S7^^.....eJ^3 = [S^^h " : (5) 

and is used to initialize a QMC simulation to calcu- 
late the effects of the Coulomb interaction U. The 



QMC result for the cluster Green function C: ^ = 
[U, G£i,...,eN thus depends on the particular disorder 
configuration {e^}. The disorder averaged cluster Green 
function is then obtained from the individual results by 

Gij — (G£j,...,ejv,!j) ' (6) 

where the average < ... >= J YliLi'i^iP{^i){- ■ ■) is to 
be taken for a system of Nc sites. Note that in principal 
QMC calculations for all possible disorder configurations 
{ei} have to be carried through. However, contributions 
from configurations with m Zn impurities on the cluster 
(m < Nc) are weighted by a factor of x™(l — x)^^"""^") 
to the integral Eq.(||). Therefore it seems reasonable for 
small concentrations x 1 and cluster sizes Nc to con- 
sider only those configurations with none or only a single 
Zn impurity on the cluster and to neglect configurations 
with more than one Zn impurity. Then the disorder aver- 
aged cluster Green function G^^ is obtained as a weighted 
sum from the two configurations as 

G^^. = xNcGl,^ + (1 - xNc)G-o,^J > (7) 

where G^^ ^- denotes the QMC result for the cluster 
Green function for the configuration with one or zero Zn 
impurities on the cluster, respectively. The prefactors fol- 
low with (1 - x)^" « 1 - xNc, xNc[l - x)^=-i « xNc for 
a; << 1 in linear approximation in the impurity concen- 
tration X. The disorder averaged cluster Green function 
Gij is then transformed back to cluster reciprocal space 
to calculate a new estimate of the cluster self-energy 

S(K) = g(K)-i -G^(K)-i, (8) 

which is then used to repeat the steps starting from 
Eq.(|) until convergence is reached. 

Results We study the superconducting instability and 
magnetic properties of the Hamiltonian Eq. (|l|) in the in- 
termediate couphng regime at 5% hole doping [n = 0.95) 
for different Zn impurity concentrations x. We perform 
calculations for Nc = 4, the smallest cluster size that 
allows for a transition to a dx2_y2-wave superconduct- 
ing state while preserving the lattice translational and 
point group symmetries. We set the hopping integral 
t = 0.25eV and the Coulomb repulsion U = W = 2eV, 
where W — 8t is the bandwidth of the non-interacting 
system. We choose V = 5eV > W + U to simulate the 
nonmagnetic closed shell ((i^°) configuration of the im- 
purity site. 

We searched for superconductivity with s-wave, ex- 
tended s-wave, p- and d-wave order parameters. As 
we found for the clean system (2^=0) |l6j, only the 
dj.2_y2--wa,Ye pair field susceptibility for zero center of 
mass momentum diverges at low temperatures. Fig.|l| 
illustrates the inverse pair field susceptibility P^^ as a 
function of temperature T for the impurity concentra- 
tions X — (circles), x = 0.05 (squares) and x = 0.10 
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FIG. 1: The inverse pair-field susceptibility versus tempera- 
ture when Nc=4,U = W^ 2eV, V = 5eV at 5% doping for 
impurity concentrations a; = (circles), x = 0.05 (squares) 
and a; = 0.10 (diamonds). The solid lines represent fits to the 
function oc (T — Tc)^. Inset: Critical temperature Tc as 
a function of impurity concentration x. 



FIG. 2: The bulk magnetic susceptibility versus temperature 
when Nc = 4, U = W = 2eV, V = 5eV at 5% doping for 
impurity concentrations x = (circles), x — 0.05 (squares) 
and X — 0.10 (diamonds). The solid lines are guides to the 
eye. The dashed line represents a fit to the function x{T) = 
Xo-XiT + C/{T-&). 



(diamonds). The corresponding critical temperature Tc 
is then calculated by extrapolating P^^(T) to zero us- 
ing the function P7^ oc (T — Tc)'*' (represented by the 
solid lines in Fig.yJ). Tc is rapidly suppressed, and the 
critical exponent 7 increases, with increasing Zn content 
X. Consistent with experiments Q this falloff of Tc is 
linear in x resulting in a critical impurity concentration 
Xc ~ 0.10, beyond which the instability to the super- 
conducting phase disappears. Compared to experiments 
{xc ^ 0.05) our calculation predicts a result about twice 
as high for the critical doping. However, considering the 
simplified description of the problem using a 2D Hubbard 
model with the complexity of the Zn impurities reduced 
to a model of diagonal disorder, we expect agreement 
only on a qualitative level. 

In order to study the correlation between the suppres- 
sion of superconductivity and possible moment formation 
by Zn substitution we study the bulk (q = (0,0)) and 
the antiferromagnetic (q = (tt, tt)) spin susceptibilities. 
Fig.^ shows the bulk magnetic susceptibility x(r) as a 
function of temperature for the same impurity concen- 
trations X used in Fig|l[ As we demonstrated in [ p^ , 
the clean system exhibits an anomaly in x(T) at a char- 
acteristic temperature T*, below which spin excitations 
are suppressed and x(T) falls with decreasing tempera- 
ture. T* thus marks the onset of singlet formation due to 
strong short-ranged antiferromagnetic spin correlations. 
In agreement with experiments [Q, the substitution of 
impurities leads to a dramatic change of the low temper- 
ature behavior of the magnetic spin susceptibility. De- 
spite the nonmagnetic nature of the impurity scattering 
term in our simulation, x{T) displays an upturn at low 
temperatures as x is increased. Moreover at x = 0.10, 



x{T) shows Curie- Weiss like behavior; i.e. it can be fit 
by the function xiT) = Xo - XiT + C/{T - 6), with 
O « 0. This clearly indicates the formation of free mag- 
netic moments due to the breaking of singlet bonds by 
the substitution of nonmagnetic impurities. 

The inverse antiferromagnetic spin susceptibility 
XApiT) for the same parameter set as in Fig.|^ is illus- 
trated in Fig.||. For the clean system (x=0) 5% hole 
doping marks the critical doping, beyond which no tran- 
sition to an antiferromagnetic state can be found . As 
can be seen from the extrapolation of the circles in Fig.^ 
the corresponding Neel temperature Tn = 0. As the 
concentration of impurities is increased we notice that 
for T ^ 0.05 XAF decreases with x. However at lower 
temperatures (T ^ 0.05) Xaf ^'^^ finite impurity concen- 
tration X falls below the curve of the clean system and 
eventually goes through zero indicating a transition to 
the antiferromagnetic state. Therefore it can be inferred 
that a finite concentration of impurities, i.e. spin vacan- 
cies, acts to enhance antiferromagnetic spin correlations. 
Our results for the site and time averaged local magnetic 
moment Txa shown in the inset to Figjs] support this 
conjecture: At low temperatures this quantity increases 
with impurity concentration x indicating that local mo- 
ments get stabilized. It is important to note that recent 
NMR measurements on YBCO [0 show a similar en- 
hancement of antiferromagnetic spin correlations around 
Zn impurities at low temperatures. 

Similar results were also obtained in |18, Eoj where an 



antiferromagnetic Heisenberg model with a spin vacancy 
was studied. In the latter report, the enhancement of an- 
tiferromagnetic spin correlations was ascribed to pruning 
of singlets by the nonmagnetic impurity in an resonating- 
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FIG. 3: The inverse antiferromagnetic (q — (7r,7r)) spin sus- 
ceptibility versus temperature when Nc = 4:, U = W = 2eV, 
V = 5eV at 5% doping for impurity concentrations x = 
(circles), x = 0.05 (squares) and x = 0.10 (diamonds). The 
solid lines represent fits to the function Xaf oc (T — Tn)^. 
Inset: The site and time averaged magnetic moment versus 
temperature; spin vacancies reduce fluctuations. The solid 
lines are guides to the eye. 



tions that mediate pairing get suppressed. Consistent 
with experiments, Tc decreases linearly with impurity 
concentration. With increasing Zn content we further 
find a change of the magnetic susceptibility to Curie- 
Weiss like behavior indicative of the existence of free 
magnetic moments. 

Our results can be understood within the RVB picture. 
The nonmagnetic impurities break local singlets and thus 
generate unpaired spins. A concomitant suppression of 
antiferromagnetic spin fluctuations that mediate pairing 
results from the pruning of RVB states Conse- 
quently, superconductivity gets strongly suppressed with 
increasing Zn content. 
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valence-bond picture (RVB). In the RVB picture, the 
short-ranged order in the doped model is described as a 
collection of nearest-neighbor singlet bonds which fluctu- 
ate as a function of time and space. When a Zn impurity 
is introduced into the system, the singlet fluctuations on 
adjacent sites are suppressed since these sites now have 
one less neighbor to form short-ranged bonds with. The 
suppression of singlet fluctuations on the sites adjacent 
to the Zn dopant also enhances moment formation on 
these sites. 

We find further evidence of the suppression of spin fluc- 
tuations by Zn impurities in the decrease of the anti- 
ferromagnetic exponent 7 towards the mean-field result 
of one with increasing impurity content x. Concomi- 
tant with this behavior the corresponding exponent 7 for 
the pair-field susceptibility (cf. Figj^) increases towards 
one. From this observation we conclude that the suppres- 
sion of superconductivity in Zn doped systems originates 
in the suppression of antiferromagnetic spin fluctuations 
that mediate pairing. 

Summary In this letter we have used the 2D Hubbard 
model with binary diagonal disorder to study the effects 
of nonmagnetic (Zn) impurities on high-temperature su- 
perconductivity. We find that Zn impurities strongly sup- 
press superconductivity. As evidenced by the increasing 
mean-field character of the pairfield and antiferromag- 
netic susceptibilities with Zn substitution, spin fluctua- 
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